On 



Exponentiation of Multiparticle Amplitudes in Scalar Theories 

M. V. Libanov, V. A. Rubakov, D. T. Son and S. V. Troitsky 

Institute for Nuclear Research of the Russian Academy of Sciences, 

60th October Anniversary Prospect la, Moscow 117312 Russia. 

(July 1994) 

Abstract 



\£> ' It is argued that the amplitudes of the production of n soft scalar particles 

(N' 

by one or a few energetic ones in theories like \(p 4 has the exponential form, 

A n ex vra! exp[-^F(An, e)], in the regime A —>■ 0, An = fixed, e = fixed, where 

00 

e is the typical kinetic energy of outgoing particles. Existing results support 
O ' 

this conjecture. Several new analytical and numerical results in favor of the 

OS 

exponential behavior of multiparticle amplitudes are presented. 
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I. INTRODUCTION 

Two classes of processes in weakly coupled theories with multiparticle final states attract 
considerable interest. One of them is the instanton-like transitions induced by collisions of 
highly energetic particles (which violate (B+L) in the standard electroweak theory) [|l],0|; 
there the typical number of final particles is expected to be of order 1/oc at energies E oc 
m/a, where m and a are mass of a relevant particle and coupling constant, respectively 
(m = m w ,a = aw in the electroweak theory; for reviews see refs. 0,|J). The other is 
the "perturbative" creation of n oc 1/A soft particles by one or a few energetic (virtual 
or real) ones in scalar theories like \{p 4 ||[5],|7|. In both cases the amplitudes calculated 
to the leading order of the perturbation theory about either the instanton or vacuum are 
unacceptably large, corrections to them exceed the leading order results and the conventional 
perturbation theory breaks down. 

In the case of the instanton-like processes, there exist strong arguments showing that 
the total probability has the exponential form [JS[-|TTH 

a oc e^ E ' E ^ (1.1) 



where Eq oc — and the pre-exponential factor is at most a power of a at E ~ Eq. Eq. fll.lD 
suggests that there may exist a semiclassical-type technique for calculating the exponent 
F(E/Eq); several ideas have been put forward in this direction |T^-[T7|. 

In this paper we conjecture that the behavior similar to eq.( |l.l|) is characteristic to 
multiparticle amplitudes in scalar theories. Namely, consider a theory of one scalar field 
with the lagrangian 

L=\{d^f-^-\^ (1.2) 

(we will set the mass of the boson equal to one wherever possible, the powers of m can be 
reconstructed on dimensional grounds). The processes we are interested in are creation of n 
scalar quanta with typical kinetic energy e by one or a few initial (virtual or real) particles, 
as illustrated in fig.0. We conjecture that in the limit 



A — > 0, \n = fixed, e = fixed (1.3) 

the amplitudes of these processes have the form 

A n oc Vri.e* F{Xn ' e) (1.4) 

where the pre-exponential factor is at most some power of A. This behavior is expected 
to be inherent in both unbroken theory (m 2 > 0) and broken theory (m 2 < 0), as well as 
in extensions of the model ( |L2|) , the function F being model dependent. If there are other 
parameters in a model, like the number of scalar fields in O(N) theory, they are all assumed 



to be fixed in the limit ( |1.3|) 



At the moment, we are unable to prove eq. fll.4p . Instead, in this paper we present a 
series of consistency checks which will be summarized in the rest of this section. Before 
doing so, let us make a few comments concerning eq.( |1.4| ). 



First, eq.( |1.4p implies that the cross section of the creation of n soft scalar particles, 

L4J2 

a n oc — ^— x (phase space) 
n\ 

is exponential in the regime ( |1.3[ ), 

a n oc e^ G(nA - e) (1.5) 

Eq. fll.5Q is completely analogous to eq. fll.lD characteristic to the instanton-like transitions. 
Whether the cross section (|1.5|) is small at all n or becomes large at n oc 1/A is determined 
by the actual behavior of the function F(n\, e) on which we have almost nothing to say in 
this paper. 

Second, there exist arguments based on the ordinary perturbation theory at low momenta 
(say, for the scalar propagator) and unitarity, which favor the exponential suppression of 
the few —>■ n cross sections at n oc 1/A ||18| . Eqs. fll.4| ) and ( |1.5| ) are consistent with this 
expectation provided that F(nX, e) is always negative and decreases with e rapidly enough. 
We think, however, that the actual computation of F(n\, e) is necessary to settle the issue. 



Finally, the obvious generalization of eq. (|1.4j) to the processes n\ — > n 2 is that the 
connected amplitude has the following form, 



A ni ^ n2 ex ^ ni \^n 2 \^ F ^ n ^ (1.6) 

which we conjecture to hold in the regime 

A — > 0, An! = fixed, \n 2 = fixed, e = fixed 



Eq.( |1.6| ) has a direct analog in the instanton-like case fL2[-f[l|; in this paper we will present 
the arguments in favor of eq. Ql.6D in the scalar theory. 

Let us now summarize the consistency checks for the conjecture ( |1.4| ). 

i) Tree amplitudes at threshold in Lp 4 theory. At e = the amplitudes 1 — > n have 
been calculated at the tree level in refs. [[^,|T^,f2"0|j , 

n— 1 

A^ ee =nl[ — \ (unbroken theory) 



/ \ n ~ 1 
A^ ee — n! ( - J (broken theory) 

In the limit ( |1.3| ) these amplitudes indeed have the form ( |1.4j ) with 

F tree = (An) b An _ An ( unbroken theory) (1.7) 

F tree = (An) b An _ An (broken theory) (1 g) 

ii) Tree amplitudes at threshold in an extension of y? 4 theory. The model of two 
scalar fields tp\ and tp 2 with unequal masses and 0(2) symmetric interaction, described by 
the lagrangian 

l = \{d^f + \{d^Y - ^ - ^jM - \^l + vl) 2 (1-9) 

has been studied in ref. | ZIfl . The tree amplitudes at the threshold for the processes 1 — > 
(ni,n2) have been calculated exactly. At n\ oc 1/A, n 2 oc 1/A they have the form 



< e l n2 oc y/n 1 \y/m\eS F ^<' niX > naX > (1.10) 

The exponent is particularly simple when ni/n 2 = 1 + 0(1/ X) where (we assume for defi- 
niteness that the incoming virtual boson is of the first type, <pi) 

F(Xn) — nX In ^ ^- 1 

V ; V 8( mi + m 2 ) I 



(unbroken theory), 



F(Xn) — nX In ^ : 1 

V ; V 2(m 1 + 2m 2 ) 



(broken theory, m\ < or m\ < 0). 

Eq. ( |1.10| ) is an obvious generalization of eq. ( |1.4j ) for the theory ( |1.9| ) . 

Hi) Energy dependence of the tree amplitudes near threshold in unbroken (p 4 
theory. At the tree level, the dependence of the 1 — > n amplitudes on A is known explicitly, 

< ree oc At (1.11) 



At small e, eqs. ( |1.4| ) and ( |1 . 1 1| ) are consistent with each other only if 



< ree (e) = A£ ee (e = 0)e An[e+o(e2)] 

where A is some constant. Since ne = E, the total cm. kinetic energy, we expect the tree 
amplitudes of the creation of n non-relativistic bosons to have the following form, 

A^ ee (E) = A^ ee (E = 0)e {AE+ - ) (1.12) 

where dots denote the terms suppressed by E/n. We calculate the leading energy correction 
in <£> 4 -theory in seci.0 of this paper and show that it indeed has the form ( |1.12|) with 

A = - 6 - 

6 

iv) Exponentiation of the loop corrections at threshold in the leading order in 

(nX). The one loop correction to the 1 — > n amplitude at threshold in (p 4 theory has been 
calculated in refs. J2"2H24f. The result is 



A l-loop = ^toce x (^ BXn 2 + (Xn)) (1.13) 

where the (complex) constant B depends on the number of dimensions. In (3+1) dimensions 
one has 

B = -— — -(In (7 + 4\/3) - in) (unbroken theory) (1.14) 

/3 

B = (broken theory) (1-15) 

At higher loops one expects higher powers of n along with higher powers of A. We show in 
sect^Q of this paper that in the fc-th loop, the leading-n correction has the form 

AT-^-(\n*) k (l + 0(l/n)) (1.16) 



with exactly the same constant B as in eqs. ( 1.14Q , ( [1.15 ). This means that these leading 



order corrections, being summed over all loops, exponentiate, 

A n = A^ ee e BXn2 x • ■ ■ (1.17) 

which is consistent with eq. (|1.4j ) with 

F{\n) = F trcc (An) + B ■ (An) 2 + • • • (1.18) 

The non-leading terms in eqs. (|1.13|) and (|1.16|) correspond, presumably, to terms of order 
(An) 3 , (An) 4 , etc., in F(Xn) and to the contributions to the pre-exponential factor in eq. (|1.4| ). 

In analogy to the instanton-like processes, eqs. (p..7|), ( |1.8|) and (|1.18|) imply that the 
tree result is actually a good approximation to F at small An, in spite of the fact that the 
one loop correction ( |1.13| ) is greater than the tree amplitude (the point is of course that 
even small correction to F leads to the large factor in the amplitude, see eq. ( |1.17| )). On 
the other hand, at An ~ 1 the series ( |1.18| ) blows up at An ~ 1, and truly non-perturbative 
techniques is needed for calculating F(An). 

v) Amplitudes at threshold in O(N) theory at large N. The amplitudes of 1 — > n 
processes at threshold in the O(N) symmetric theory of N scalar fields have been calculated 



in ref. [25] to all loops in the regime 



A — > 0, XN = fixed, n = fixed 



;i.i9) 



They are equal to 



-4 



A' 



A 



n/2 



n i 



;i.2o) 



\8(l + cXN)J 

where the coefficient c depends on the renormalization scheme. The regime ( |1.19[ ) is different 
from the regime 

A — > 0, Xn = fixed 

with all other parameters (including N) fixed, which is considered throughout this paper. 
The two regimes match only in the region nX <C 1,NX <C 1 where eq.( |1.20| ) is consistent 
with eq.( p..4| ) with A^-dependent pre-exponential factor 

aN const- (An) N _ Atree,N=l 



So, there is no contradiction between the large-A^ results and our conjecture. 

vi) Infrared behavior of the amplitudes in (p 4 theory in (2 + 1) dimensions near 
threshold. Loop contributions to the amplitudes exactly at the threshold are infrared 
divergent in the (p 4 theory in (2 + 1) dimensions. At small but non-zero e, the terms of order 
(Alne) fc appear in fc-th loop. These terms can be summed up by the renormalization group 
technique p6| with the result 



n{n — 1) 



A.M=4r(i-±m e y = 



(unbroken theory) 



;i.2i) 



x(n-l) 



A n (e)=AT e (l + ^lacj 
which is valid in the leading logarithmic regime 



(broken theory) 



A — > 0, e —> 0, A In e = fixed, n = fixed 



(1.22) 



This regime matches with one considered in this paper, eq.( |1.3|) , at 



e< 1, Alne< 1, An < 1 (1.23) 



where eqs. ( |1.21| ) and ( |1.22j ) agree with eq.(|1.4|) with 



3 

F(e) = F trce H (An) 2 lne (unbroken theory) 

327T 

F(e) = F tiee - —(An) 2 lne (broken theory) 

Similar results hold for O(N) extension of the y? 4 theory. In the unbroken case one finds in 



the limit ( |L23D and fixed N 



F(6)=F^ e + ^£hir 



itrcc _ 

327T 



We present the leading logarithm calculations of the 1 — > n amplitudes in (2+1) dimensions 
in sec£.[7T/[ of this paper for completeness. 

vii) Finally, the conjecture ( |1.6[ ) on the behavior of the connected amplitudes rii — > n 2 
can be tested at the tree level. The dependence of the amplitudes on A follows from counting 
the vertices in the tree graphs, 

A%% n2 oc X^ 



so one expects from eq. ( |1.6|) the following functional form of the tree amplitudes at rest at 
large n x and n 2 , 

A ni ^ n2 = n 1 !n 2 !A Iil ^e (ni+n2),J>(n2/ni) (1.24) 

where $(n2/ni) is some unknown function. We present numerical results confirming the 
scaling behavior, eq.( |1.24|) , in sect!^ of this paper, where heuristic arguments in favor of the 



conjecture ( |1.6| ) are also given. 

The rest of this paper is devoted to the evaluation of the new results entering the above 
arguments: energy dependence of the tree I —> n amplitudes near the threshold (sect.y), 
leading-n corrections to 1 — > n to all loops at the threshold (sect.|lll[), leading-log infrared 
behavior of 1 — > n amplitudes to all loops in (2+1) dimensions (sect.[IVD, and n\ — > n 2 tree 

8 



amplitudes at rest at large ni, n 2 and fixed r^/ni (sect.0). Although these new results 
may be of interest by themselves, we think their main value is to support our conjectures, 
eqs.flOJ) and (pp. 

Sect.[VI| contains concluding remarks. 

II. ENERGY DEPENDENCE OF MULTIPARTICLE TREE AMPLITUDES 

AROUND THRESHOLD 

Our first consistency check of the exponentiation hypothesis is the calculation of the 
multiparticle tree amplitudes when final particles are not at rest, but have finite kinetic 
energies. We will be interested in the amplitude of the production of n final particles with 
momenta Pi from an initial virtual particle in the y? 4 theory without symmetry breaking. In 
what follows we consider the case of non-relativistic final particles, Pi <C 1. Denote the total 
kinetic energy of final particles in the center-of-mass frame by E, in this frame (J2p% = 0) 
one has E — ^J2pf ( we se t the mass of scalar bosons equal to 1). More generally, one has 
the following formula for E in arbitrary frame, 

in i n 

in in 

1=1 l^J 

The case of threshold production corresponds to E = 0. In this section we show that at 
large n and E -C n, the tree amplitude has the following form 

71-1 

A n (j>i,...,p n )=n\(~\ 2 e~l E (2.2) 

As the first step for establishing this formula, we consider the region of very small E, E <C 1, 
where the amplitude can be calculated for arbitrary n, not necessarily large. 



A. Amplitude at E <C 1 and arbitrary n. 

If the momenta of final particles Pi are small, we can expand the tree amplitude A n (pi) 
in powers of Pi. The zeroth term is the threshold amplitude A n . Since the amplitude is P- 
invariant, A(pi) = A(—pi), the first energy correction to the threshold amplitude is quadratic 
in pi. The Galilean invariance and the symmetry of the amplitude under permutations of Pi 
ensures that this quadratic correction is proportional to E. So, one writes 

A n (p h ..., Pn ) =A n + a n E + 0(E 2 ) 

where a n is some constant that depends on n. The purpose of this subsection is to calculate 

OL n . 

The tree amplitude A(pi, . . . ,p n ) satisfies the following recurrence relation (see fig.|2|, cf. 

((n + E) 2 - (f» 2 - l)A n ( Pl , ..., Pn ) = 
i=i 

A E E-V(pf ) ,... ) pg)A tt (pf ) ,...,pg)^.(pS 8) ,...,p£ ) ) (2-3) 

ni,n 2 ,n 3 -p 

where the sum is taken over all n\, n 2 , n 3 satisfying the relation n\ + n 2 + n% = n and over 
all permutations V of momenta pi. In fact, the sum runs over odd values of n\, n 2l n^ only, 
since in the theory without symmetry breaking the 1 — > n amplitude vanishes at even n. 
For convenience we will work in the center-of-mass frame where Y%=iPi = 0- Hereafter we 
include the initial propagator into the amplitude A n (pi, . . . ,p n ). 

Let us expand both left and right hand sides of eq.( [2.3|) in powers of E. At the zeroth 
order in energy, one obtains the recurrence relation for the threshold amplitude A n , 

TV 

(n 2 -l)A n = X E , '-. , A ni A m A n3 (2.4) 

while in the first order in energy E one has 

Tl Tli — 1 

(n 2 - l)a n + 2nA n = 3A E . ', . -, a ni A n2 A n3 (2.5) 

n ll n 2 ^;n 1 >l"l !ri 2!n 3 ! 71-1 
10 



By introducing the generating functions 



the recurrence relations, eqs. 
for A(t) and a(r), 



oo A 

Mr) = £ %e nT 



n=l 



n\ 



oc 






n=3 



(2.6) 



J|), can be rewritten in the form of differential equations 



{d 2 - 1)A(t) = XA 3 {t) 



{d 2 T - l)a(r) + 2d 2 A(r) - 2d T A(r) = 3Xa(r)A 2 (r) 



(2.7) 



In particular, the equation for A(r) coincides with the field equation for spatially homoge- 
neous configurations. Let us discuss the boundary conditions for A{r) and a(r). The 1 — ► 1 
amplitude at arbitrary momentum is A\ = 1. So, the boundary conditions are 



A{t) = e r + 0(e : 



3t\ 



q (t) = 0(e 6T ) 
in the limit r — > — oo. The solution to eqs.( [2.7| ) with boundary conditions 

A{T) = T^ 



Oi(T 



5A 



,3r 



1- ^e 



A„2t 



(l - I e- 



Expanding ( [2.9|) in powers of e r , we obtain A n and ct r , 



K - »! m ' 



n odd 



is 



(2.8) 



(2.9) 



«, 



_„,^V i( B __L ,, 

1 8/ 6 V n-1," 



n odd, n > 3 



11 



Thus, up to the first order in E, the tree amplitude is 



A(pi,...,p n ) = n\ 



n— 1 



E 



(2.10) 



\6 6(n-l) 

From eq. (|2.10| ) one sees that the expansion in E is a good approximation only when S<1, 
i.e when the total kinetic energy of final particles is much smaller than the energy required 
for creating an additional scalar boson, and the energy per a final particle is much smaller 
than 1/n. However, at large n, one can calculate the amplitude in a larger range of E. In 
the next subsection we will obtain the expression for the tree amplitude at E <C n, which 
covers the whole non-relativistic final phase space. 



B. Non relativistic tree amplitudes at large n 

For reasons explained in Introduction, we expect that at large n and E -C n the amplitude 
has the following form, 



A(p 1 ,...,p n ) = A n e 



■»* = n![^l " e AE 



(2.H) 



where A is some constant. Let us verify that at large n the Ansatz (j2.ll ) satisfies the 
recurrence relation (|2.3|) with the accuracy of E/n <C 1. 

Substituting (|2.11|) to (|2.3| ), one has on the left hand side of eq.(|2] 



((„ + Ef - l)A(pi) = n 2 A n e Ah (l + 0(E/n)) 



(2.12) 



Consider now the right hand side of eq. (|2.3|) . Taking into account eq. (|2.1| ) , one writes 



A 2_^ 7 t A ni A n2 A n:i exp 

ni,ri2,ri3 -p 






i¥=j 



U: v ^ n.i!n 2 !n 3 ! ^^_ _ ( A ^ (1) 2 A 



SA n e AE £ 



711,712,713 



n! 



53 exp 



7^ 



2ni 7T-! J 



(2), (3)) 



(2.13) 



Let us expand the exponential function in powers of momenta. The first term is equal to 
one. If one retains only this term, the sum is 



12 



ni\n 2 \rfi\ n 2 - 1 



77 ' 
ni,ri2,n3 ~p 



(note that the sum runs over odd values of rii, n 2 , n 3 ). One finds that the expressions fl2.12| ) 



and fl2.13| ) coincide with the accuracy of 1/77, i.e., the recurrence relation ( |2.3| ) is satisfied 



with this accuracy. So, to show that eq. (|2.3| ) is satisfied by our Ansatz, one has to check 
that the contributions from higher terms in the Taylor expansion of the exponent to the 
sum in eq. (|2.13| ) are suppressed. 



Let us check that the second term in the Taylor series of the exponent indeed gives a 
small contribution to the right hand side of eq.( [2.3| ). To evaluate the sum that comes from 
this term, 

E ^ E (£ Erf' 2 + £ Ep! 1 >«-?> + (d) - (2), (3)) 

one makes use of the following combinatoric relations, 

ni!n 2 !n 3 ! ^^ ( 1} 2 774 " 2 

"'■ V i=\ i=\ 

77i!77 2 !r7 3 ! ^^ (i) (i) ni(ni-l)A 
i 2^ 1^ Pi Pj = — ( TV 1^ P'Pi 



The contribution from the second term in the Taylor series to (|2.13|) is then 



-8A n e AE V -W+ ? ^ ifp ift 
™V2r7^ 77(77-1) Jf^ 

= -16A„e A£; (77 + l)AE (2.14) 

where we have made use of the relations J2p 2 — 2-E, J^PiPj — ~E. It is clear that at large 
77, the expression ( [2.14|) is suppressed by a factor of E/n as compared to the leading term, 
eq.( |2.12| ). So, in the case E <C n the recurrence relation is satisfied with accuracy of 0(E/n) 



by the ansatz ( p. 11 



The recurrence relation at large 77 does not determine the constant A. To obtain the 
value of A, one can consider the region of very small E, E <C 1 and make contact with our 
previous result fl2.10|) which in the large n limit reads 



13 



'A\ ( - 1)/2 / 5 



A n (E)=n\[-\ \ 1 -q E ) ( 2 - 15 ) 



By comparing (|2.11|) and ( |2.15|) , one finds 

A=- S - 

6 

So, we have established the formula fl2.2|) for the tree amplitude of the production of n non- 
relativistic particles. The result has the form expected from the exponentiation hypothesis. 

III. EXPONENTIATION OF LEADING -n CORRECTIONS FROM ALL LOOPS 

FOR I^tjAT THRESHOLD 

A. General formalism 

In this section we consider the process 1 — > n (n odd) at threshold in the unbroken tp A 
theory. The amplitude of this process has been calculated in refs. [0,0,^] at the tree level 
and in refs. [p2|-|2"4]| in one loop. The result is 

Al ee + A^ loop = A^ ee (l + B\n 2 ) 

where B is a complex number that depends on the number of space-time dimensions. In 
(3+1) dimensions the value of B is given by eq.( [1.14| ). We will consider higher loop correc- 
tions to the leading order in n. In the fc-th loop they are of order (Xn 2 ) k ||24|| . We will show 
that these corrections sum up into the exponent 

A n = Al ee exp (BXn 2 ) 

This result strongly supports the conjecture stated in Introduction. 



The technique that we use in this section is based on the formalism of ref . [Bj that reduces 
the problem of the evaluation of the amplitude 1 — > n at the threshold to the calculation of 
one-leg Feynman graphs in a certain spatially homogeneous classical background field. We 
perform direct evaluation of Feynman graphs, and show that to the leading order in n, one 

14 



can reduce the calculation of loop graphs to a simpler calculation of tree graphs in some 
effective theory. 



Let us outline briefly the technique developed in ref. [ I5|| . For convenience let us again 
set the mass of the particle m — 1. Consider the transition from an initial virtual particle 
with four-momentum P M = (n, 0) into n final particles, each with four-momentum (1,0). 
The reduction formula for the amplitude can be written in the following form, 

d n r 
A n = Jim limy - l) n g^J rf 4 xe-^(0|^(x)|0) J=JoexpM) , (3.1) 

where the matrix element is calculated in the presence of the source J = Joexp (iut). 
Let us write the classical field equation with the source, 

dfa-ip- \<p 3 + J Q e wt = 0. (3.2) 

One of the solutions to this equation has the following expansion, 

<p (t)=z{t) + ..., (3.3) 

where 

z(t) = z e u = -^-e u , (3.4) 



and dots stand for terms proportional to e 3lt , e 5lt , etc. The crucial point is that one can 
take limits uo 2 — > 1 and J — ► simultaneously, so that the ratio Zq remains finite. In this 
case, the field equation Q3.2|) becomes sourceless and its solution having the form (|3T3| ) can 
be found exactly, 



*(*) 
1 - \z\t) 



Mt) = ; I' (3.5) 



Taking into acount the following relation which comes from eq.( |3.4| ), 

one can rewrite the amplitude in the following form 

15 



^ = ^<%l°> (3- 6 ) 

where the matrix element is calculated in the classical background ( |3.5|) . 

Thus, the problem of evaluating the 1 — > n amplitude at the threshold reduces to the 
calculation of the matrix element of the field operator in a classical background. This matrix 
element can be calculated in terms of conventional Feynman graphs. This is the technique 



invented in ref. |19 



At the tree level, one replaces the matrix element in eq.( |3.6| ) by the background field 
( p.5| ) and finds 



AT e = n! I - I " 



n 



To calculate loop corrections, it is convenient to evaluate the matrix element in euclidean 
space and after that perform the Wick rotation. By introducing the euclidean time variable, 

1. A . .71 

T — it H — m — h lnz Q + i— , 

2 8 2' 

the background configuration in the euclidean space is written as follows, 

f\ 1 



2 cosh r 
Extracting the quantum part (p from the field operator <p, 

if = tpo + (p 

one has (0|<^|0) = tpo + (0|^5|0). To calculate the matrix element of (p one first derives the 
Feynman rules. For this purpose we write the euclidean lagrangian for (p, 

from which one obtains the Feynman rules shown in fig.|3]. 

Let us denote the contribution from the k— th loop by ipk{ T )- Since the background field 
cpo ~ cosh _1 (r) is singular at r = in/2, one expects that ipk is also singular at this value 

16 



of r. Assuming that the leading singularity is a pole of order n k (in fact, we will see that 
n k = 2k + 1), one can expand <£>& around the singularity in powers of ip , 

Substituting eq.fl3.5|) into eq.fl3.7p and differentiating n times with respect to z, one finds 
that the leadings behavior of the amplitude is determined by the leading singular term in 
the right hand side of eq.( [3.7|) . Namely, at large n the contribution from the fc-th loop to 
the amplitude will be 



Co 



- 2 n (n fc -l)/2 



(-) n nk - l + 0{n nk ' 2 ). 



K-l)! \\. 

So, if we wish to calculate, in each order of the perturbation theory, the leading-n contribu- 
tion only, we do not need to know exactly the function (fj, but only its leading singularity 
at r = 27r/2. 

At the one-loop level, the only graph that makes contribution to (0|<£>|0) is shown in 
fig.f|a. Analytically, one has 

ip x (x) = (-3A) [dyD(x,x')(p (x')D(x',x'), (3.8) 



where D(x,y) is the propagator in the background <^o- Let us consider it more closely. 
D(x,x') satisfies the following differential equation, 

OD(x, x') = {-dl + 1 + 3\<f 2 )D(x, x') = 5\x - x'). 



The analytic formula for the propagator is known |22|j2J| . It is convenient to write the prop- 
agator in the mixed coordinate-momentum representation (momentum in space, coordinate 
in time), 

d d p 



D(x ' X,) = /(2^ eiP(X_X ' )jD(r ' T/;p) 



(2tt) c 

where d is dimension of space and r and r' are time components of x and x f , respectively. 
One writes 
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D(t, t>; p) = ^r(f?(T)fZ(T')e(T> - t) + mr)f?(T')9(T - r')) (3.9) 

where ff{r) and /^(r) are two linear independent solutions to the homogeneous differential 
equation 

6 



-di + u 2 - 



cosh T 



) f(r) = 



with oj = -\/p 2 + 1- /^(t) and /^(r) tend to zero as r — ► — oo and r — > oo, respectively, 

3 



/j"(r) = e" T u 2 - 3wtanhr + 2 — 

cosh r 



f%{r) = /r(-r) 

In eq.(^) W p = /1/2 — /2/1 is their Wronskian, 

W p = 2cj(cj 2 -I)(cj 2 -4). 

In what follows it will be more convenient to work with new functions / and g which are 
linear combinations of the old ones, 

fi{r)=^l 2 {f{r)+g{r)) 

h{r) = e— / 2 (/(r) - g(r)) 

The functions / and g have different behavior around r = in/ 2: while / is singular, /(r) = 
—3 cosh - r, the function g is regular at this point, g(r) ~ cosh r. In terms of / and g, the 
propagator can be represented in the following form, 

D{t, r'; P ) = A)(t, r'; p) + D x (r, r'; p), 

where 

^o(r,r';p) = ^/(r)/(r') 
and 



Di(t,t>;p) = -U(r - r')(f(T)g(r') - f(r')g(r)) - g(r)g(r')]. 

Do contains the strongest singularity of D, while D\ is less singular than Dq. Note that 
Dq(t, t', p) factorizes, this fact will be extensively explored in what follows. We will need 
the following formula for D which is correct in the sense of the leading singularity, 

^(ryjp^Jr^TKV) (3-10) 



Let us consider the propagator at coinciding points that enters eq. (|3.8| ). The leading 
singularity in D(x,x) is 

D(») S D( X , ») 2 /|* ^-^-1^-^f (T) 

/"_^P 9 1 a2 ^ 4, x 

7 (2vr) d 2cu(cu 2 - l)(cu 2 - 4) cosh 4 (r) ^ ' 

where we introduced the notation 

f d d p 9 

~ J (2^p 8w(cj 2 -1)(^ 2 -4) ^ ^ 

Now we are able to calculate <pi. Acting by the operator O on both sides of eq. (|3~lf ) and 
taking into accound the relation OD(x,x') = 5(x — x'), one finds 

6(p! = -3\ 3 B(p 5 (3.12) 

To solve this equation let us notice that 

o^ (si + 1 - -4^ wr /2 ' 



cosh tJ \XJ cosh™ 



2\ n / 2 (n-2)(n + 3) 

Xj COsh n+2 T 



r( T ) ^— r^-^ + 0(1/ cosh" r) 



= — I — 1 (n — 2)(n + 3)v9q +2 + subleading terms 

where only the term with the strongest singularity is written explicitly. Reversing this 
equation we obtain 
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(6-Vo +2 )W= d d yD(x,y)^ +2 (y) 



\o (3-13) 



(n-2)(n + 3)A 
Applying this formula to eq. fl3.12j ) one finds, in the sense of leading singularity. 



(pi = X 2 B V 3 (3.14) 

So, the amplitude up to one-loop is 

A n = A^ ee {l + BXn 2 + 0{Xn)) 



which coincides with the result obtained in refs. |22|,|24[], eq. fll.l3| ) 



So, we have seen that the one loop correction can be evaluated without cumbersome 
calculations. Let us turn to the two-loop level. 

B. Two loop level 

In the previous subsection we have calculated the tadpole of fig.|]a. In our calculation 
we have replaced the propagator at coinciding points, D(x',x'), by its leading singular part 
which is proportional to <Pq. For further convenience, we represent the amplitude after this 
replacement by the graph shown in fig.|]b, which can be obtained from fig.f|a by cutting the 
inner line. Each line that ends on a bullet corresponds to a factor of B l l 2 X(p 2 ) . 

There are six different graphs that contribute to the two-loop amplitude (see fig.|5|a-f). 
The first two graphs are easy to calculate by using the same technique as in the one-loop 
case. One writes for the first graph (hereafter equalities are to be understood in the sense 
of leading singularities) 

¥2a = -3AO~Vo<Pl 



Substituing cpi from eg. ( \i. 14 ) into this equation and applying the formula ( |3.13 ), one obtains 

1 



v4 D 2,„5 



p 2a = -X^BYo 
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Analogously, for the second graph, fig.§b, one has 

cp 2b = -QXO- 1 i Pl D = h A B 2 if 5 

Graphically, one can represent the first two graphs in the new form as shown in the first two 
equations of fig|J 

Let us consider the third graph. One writes 

(0<p 2c )(x) = 18X 2 (p (x) JdyD 2 (x,y)(po(y)<pi(y). (3.15) 

At first sight, to find the leading singularity in (p 2c , it appears to be natural to suppose that 
D(x, y) can be replaced by D (x, y). In this way one would have 

(Oip 2c )(x) ~ Mr)f (r) j -^ — 2 (3.16) 

The integral is some constant, so one would obtain 

(6<f 2c )(x) ~ <^o(^) 

where we have made use of the relation /(r) ~ V 9 o( r )- ^o, one would conclude that (p 2c ~ 
<Pq which is less singular than the result for the graphs |5]a and |5|b, (f^. However, this 
argumentation is wrong. The reason is that D q (t,t') is smooth at r = r', so that the 
contribution proportional to D$ into the right hand side of eq. fl3.15D is not too singular. 
More singular contribution comes from the term D Di, 



{6(p 2c ){x) = 36\ 2 ip (x) j dy D Q (x,y)Di(x ) y)(p G {y)ipi{y). 

Since the contribution with D 2 in the integrand is more regular, one can replace D\ by D 
and not spoil the leading singular behavior, 

(Oip 2c )(x) = 36A 2 ip (x) j dyD (x,y)D( y x,y)ip (y)<p 1 (y). (3.17) 

Switching to the mixed representation and using eq.( |3.10| ), one obtains 

(6 V2c )(x) = 3Q\V (x) l-^- 3 dr^D(r,r';p)^(r f )^(r f ) 
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Let us first take the integral over dr' . Since one is interested only in the leading singular- 
ity, one can replace D(t, r'; p) = (— d 2 + u? — 6 cosh -2 r)~ 1 (r, r') by O -1 = (— d 2 + 1 — 
6cos1i~ 2 t) -1 , since 1 and u> are both more regular than d 2 and cosh _2 r. The integral over 
d 3 p gives an additional factor (— B), so one has 

ip 2c = -36A 4 Sd"VoO~VoVi ( 3 - 18 ) 



Making use of the formula for tpx, eq. (|3.14j ), and the rule ( |3.13|) , one obtaines finally 

^ = ^B 2 ^ 



Let us now turn to the graphic interpretation of eq. (|3.18|) . The corresponding graph is 



shown in fig.|^c'. The solid line represents the operator O" 1 , which originates from the full 
propagator in fig.|5|c', while each line ending with a bullet comes from cutting a propagator 
and corresponds to a factor of \B 1 / 2 ^. The additional factor 2 comes from the two possible 
internal lines we can cut. It is instructive to learn which graph corresponds to eq.( [3.16| ). 
This graph is shown in fig|| it can be obtained from the initial graph by cutting both 
internal propagators. Fig|5] is a disconnected graph, and it is clear that it is proportional 
to the one-loop graph, fig.|]b, i.e., it is less singular than any typical two-loop contribution. 
Hence this graph can be neglected. 

Now it is clear how to evaluate the remaining graphs in fig|| One must cut maximum 
number of internal lines in such a way that the obtained graph is tree, but still remains con- 
nected. If a graph can be cut in different ways, all the obtained graphs make contributions. 
For instance, for the graph of fig.^f, there are five different ways to cut: one can cut the line 
1 and 3, 1 and 4, 2 and 3, 2 and 4, or 3 and 4. 

To make our presentation more explicit, we write down here the result for all two-loop 
graphs (fig.^D 

<P2a = -(-6A)6~Vo<^ 2 , 

<P2b = ^(6X 3 B)d-V fi, 
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<^ 2c = -36A 4 5CrVoO Wi 



cp 2d = i(36A 6 i? 2 )6-V^~V^ 



V2e = (6X 3 B)d- 1 V i ¥>i, 



<p 2f = (-6A)6-Vo^ - 2(36A 4 5)CrVoC r Vo<^i- 

Summing up the contributions from all graphs, one obtains the total two-loop correction to 
the generating function, 

^ = !(-6A)d-Vo^ + ^(6A 3 5)0-V^i - 3(36A 4 J B)6-VoO-Vo^i+ 

3^(36A 6 S 2 )6-V^-Vo- 
o 



Making use of the explicit formula for ip\ and the rule (|3.13|), one finds 



y? 2 = 3A 4 £^>. 

(We have checked this result by direct — and tedious — calculation of the Feynman graphs.) 
Thus, the leading-n threshold amplitude up to the two loop level is 

< ree + A^ oop + A 2 - loop = < ee (l + BXn 2 + ^A V j , 

We see that, up to the two-loop level, the leading-n correction coincides with the expo- 
nent exp(5An 2 ). In the next subsection we show the exponentiation in all orders of the 
perturbation theory. 

C. Exponentiation 

In the previous subsection we have designed a technique that reduces the calculation 
of the leading-n behavior of any Feynman graph to a simple arithmetic operations. The 
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prescription is to cut a maximum number of internal lines to obtain a set of connected tree 
graphs with external lines. Each external line that comes from an internal propagator is then 
replaced by a factor of XB 1 ^ 2 ^, and each propagator that has not been cut is associated 
with the operator O -1 . This procedure allows us to calculate the leading singularity in an 
arbitrary Feynman graph. 

Consider the fc-loop level. Since the new graph obtained by cutting is tree, it is clear 
that one must cut k propagators. So, the constant B appears k times. By a simple counting 
argument one can show that every fc-loop graph is proportional to \ 2k B k tp + . Thus, the 
generating function to all orders of perturbation theory can be written in the following form 

oo 

(0M0) = ^ ]T4(A 2 I^) fe (3-19) 

fc=0 

where dk are some constant. 

According to eq.( |3.6| ), this formula implies the following leading-n expansion for the 
amplitude A n , 

(00 ofc J \ 

1H-£^W*J (3.20) 

To see that the whole sum (|3.20p is in fact the exponent, we have to show that dk = 



(2k)\ / (2 k k\) . There are at least two ways to see that this is indeed the case. The first one is 
indirect and is based on the fact that the coefficients dk do not depend on the dimensionality 
of space-time (in our techniques, the dimensionality of space-time enters only through the 
constant B that is determined by eq.( |3.11| )). So, to calculate the constants dk one can work 
in any space-time dimension. One notices that in (2+1) dimensions the integral that defines 
the constant B, eq. (|3.11| ), diverges logarithmically in the infrared region (see section |TV| ). 



One can regularize the integral by introducing, for example, finite spatial volume V. In this 
case the series ( p . 1 9| ) is in fact a sum of leading logarithms. Comparing this series with the 



leading-log result of section [TV], one finds that the right hand side of eq. ( |3.20| ) in fact sums 
up to an exponent. 

The second way is a direct analysis of the perturbative series. First, we note that the 
graphs obtained after the cutting procedure have the same form as the tree graphs in an 
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effective theory where the condensate y?o is shifted by A-B 1 / 2 ^, with the only difference in 
the symmetry factor: each graph with 2k external lines ending at bullets must be multiplied 
by a factor of (2k)\/(2 k k\). So, the prescription is as follows. One searches for a solution to 
the classical field equation (which is equivalent to the summation of tree graphs), 

%<p = Xcp 3 (3.21) 

(the mass can be neglected) which has the form 

<p cl = <fo + AB^Vo + • • • 

The meaning of ... in this formula is explained in the explicit expansion of tp c i, 

oo 

<p d = <PoY, a ^B 1/2 <Po) k (3.22) 

fc=0 

where a^ are some numerical coefficients, «o = «i = 1- For k > 1 the coefficient a^ can 
be calculated perturbatively and is given by graphs with k external lines ending at bullets. 
Since dk are also given by graphs with 2k lines ending at bullets, and because of the difference 
in symmetry factors in the two theories, this analogy implies the following relation between 
the coefficients dk in eq.( |3.19| ) and a^ in eq.( |3.22j ), 

(2*0 ! 

4 = Wr a2fc - (3 - 23) 



We are interested in the region around the singularity r = wr/2, so in eq.( |3.22 
replace tp by its singular part 

'2 1 



one can 



fa 



A r-m/2 



and obtain 



V A r — lit/ A (t — m J iy jz: 2 \ \ A r — nx / 2 J 

where we have written down explicitly the first two terms of the sum. On the other hand, 
one can verify that 

25 



' e[T) = J r-in/L^XB (3 ' 25 ' 

is a solution to eq.( |3.21| ) and it has the same form as indicated in ( |3.24j ). Comparing 
eqs.(^2|) and ( |3^4|) , one obtains a k = 1. From ( gggfl we see that 4 = {2k)\/{2 k k\), so the 
leading-n amplitude, eq. fl3.20|) , is indeed exponential. 



So, we have shown that leading-n corrections sum up into the exponent. We note finally 
that the whole calculation of this section can be applied to the y? 4 model with reflection 
symmetry breaking without any major modification. The leading-n loop corrections also 
sum up into the exponent, but in the broken case the constant B is real and positive. So, 
instead of being reduced, the amplitude is enhanced by a factor of exp(B Xn 2 ) . Whether this 
enhancement persists when one includes more terms in the expansion of the function F(Xri) 
is still an open question. 

IV. RENORMALIZATION GROUP FOR MULTIPARTICLE PRODUCTION IN 
(2+1) DIMENSIONS NEAR THE THRESHOLD 

A. Unbroken <p 4 theory 

In this section we consider the multiparticle production in (2+1) dimensions, where the 
infrared divergencies exist in a certain kind of loop graphs, which break the conventional 
pertubative expansion in a region close enough to the threshold. The graphs of this kind 
contain loops related to the rescattering among final soft particles. It is a peculiar feature 
of (2+1) dimensions that right at the threshold, rescattering graphs diverge logarithmically, 
so even at small number of final particles n, the calculation of the amplitudes around the 
threshold requires a nontrivial summation of an infinite set of graphs. 

We describe a technique to perform this summation which is a modification of the conven- 
tional renormalization group (RG). This technique allows us to sum up leading logarithms 
from all orders of the perturbation theory. By this technique we obtain the amplitude 
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1 — » n in the y? 4 models with both broken and unbroken reflection symmetry, and in the 
0(N) theory 

Let us consider the (p A theory in (2+1) dimensions, 

where, as before, we set the mass of the particle equal to one. This is the theory of interacting 
relativistic bosons. 

To describe bosons in the low-energy limit, one writes the following effective lagrangian 
in terms of a non-relativistic bosonic field \l/, 

L cS = ¥id y - i(<9^t)(<M - g&tfm, (4.1) 

where g is some, yet to be determined, effective coupling. 

Non-relativistic bosons, interacting via a delta-like potential (as in eq. (|4.1| )), have been 
known for long time as an example of a non-relativistic system with dimensional transmu- 
tation |p7|] . In fact, the counting of dimensions in (2+l)d (in non-relativistic kinematics 
the relation between the dimensions of energy and momentum is [E] = [p] 2 ) implies that 
\j/t\[/tv|/v]> j s a marginal operator and the coupling constant g is dimensionless. Apparently, 
there is no scale parameter in the theory described by the lagrangian (|4.1| ). However, this 
scale exists and is merely the boson mass (which is equal to one in our notation), which 
plays the role of the "ultraviolet" cutoff in the effective theory. 

To make contact between the effective lagrangian and the initial Lorentz-invariant one, 
one compares formulas for the amplitude of elastic scattering of two bosons computed in 
both theories. This results in the following relation between g and A, 

g = 3A/8. (4.2) 

Considering g as the bare coupling, or the coupling at the scale of the ultraviolet cutoff, 
one can study the evolution of g as a function of the momentum scale. For this purpose 
we introduce the running coupling constant g(t), which has the physical meaning of the 
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strength of the interaction between bosons at the momentum scale p = e *, (t — ln(l/p)), 
or, equivalently, the kinetic energy 



e" 2t 



There is only one Feymnan diagram (shown in fig.0) that makes contribution to the corre- 
sponding beta function. Simple calculation yields the following RG equation, 

dg(t) = g 2 {t) 
dt it 

which has the solution 

/ n Q 3A / 3A \ , 

9W = TT|* = T( 1 + S i ) • (4 ' 3) 

where we have made use of the initial condition for g(t), g(0) = g, and g is defined by eq.( |4.2| ). 
From eq. ( |4.3|) one sees that the strength of the interaction between bosons decreases as the 
momenta of particles tend to zero. Later, we will demonstrate that this property holds for a 
more general case of O(N) model, while in the theory with symmetry breaking the behavior 
of g(t) is just the opposite. 

The coupling constant changes considerably from its initial value only when the momen- 
tum scale is exponentially small in A, i.e., t = ln(l/p) ~ 1/A. So, the renormalization group 
is suitable for considering the regime A — ► 0, Xt ~ 1. All further considerations in this 
section will be done for this particular regime. 

In fact, the flow of the effective coupling, eq. (|4.3| ), can be obtained by a simpler method 
of direct summation of bubble graphs: in non-relativistic theories only these graphs con- 
tribute to the elastic scattering 2 — > 2. However, in the calculation of the amplitudes of 
multiparticle production near the threshold more complicated diagrams are involved and 
the problem cannot be reduced to the summation of bubble graphs. For example, for the 
1 — > 3 process the one- and two-loop diagrams that make contribution in our leading-log 
regime are presented in fig.[8]. To deal with these processes the RG technique is essential. 
In the diagrammatic language, the renormalization group corresponds to the summation of 
leading logarithms, i.e., terms proportional to (Xt) k in the whole perturbation series. 
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Let us try to describe the production of n final particles in terms of the non-relativistic 
bosonic creation operator W . The only relevant candidate is A n fy^ n . The following relation 
should take place, 



7? 



<^|0> = (n\A n ^ n \0) eS , 



where the matrix element in the left hand side is written in the initial <£> 4 model and the 
right hand side is understood as a matrix element in the effective non-relativistic theory. 
From this equation one finds that A n is equal to the 1 — > n amplitude when the spatial 
momenta of final particles are small, but not exponentially small (so the logarithms do not 
appear in loops). Let us for simplicity consider small enough n where A n coincides with the 
tree 1 — » n amplitude |7| jr9|j20|] , 

/ A \ (n-l)/2 

A n = AT C = n! M • (4.4) 

When the characteristic scale of momenta of final particles is exponential, the right hand 
side of eq. (|4.4|) is substantially renormalized by loops. One can treat this renormalization 
by introducing the running vertex A n (t) and solving the corresponding RG equation, 

dA n (t) n{n - 1) 1 



dt 



IX 



: g(t)A n (t). (4.5) 



Graphically, this equation is represented in fig.|9|. Substituting the function g{t) given by 
eq. (|4.3|) , to eq. (|4.5|) , we obtain the function A n (t), 



w(n-l) n(n-l) 

A n (t)=A n (l+ 9 -f-t\ 2 =< rec (l + ^J 2 • (4-6) 

So, we have found that the RG technique allows us to calculate the 1 — > n amplitude A n (t) 
in a region close to the threshold where Xt ~ 1. Note that exactly at the threshold, i.e., at 
t = +oo, the amplitude vanishes. 

B. Broken </? 4 theory 

The above analysis is equally applicable to the case of broken symmetry, 
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with the only exception that for evaluating the bare coupling g entering into the lagrangian 
of the effective theory one should take into account not only the diagram with a four-boson 
vertex, but also diagrams with two three-boson vertices, i.e. those with the exchange of a 
virtual particle in s— , t- and u-channels. In contrast to the unbroken case, the resulting 
amplitude is negative, which means the attractive character of the force between bosons at 
low energies. One has 

3A 
9 = ~Y 

(we assume v 2 = (2A) _1 , so the mass of the boson is equal to one). Eq.(|4.3|) implies that 

the strength of the interaction increases with t, 

Q 3A / 3A \ ., . 

»W = TT|7 = "T I 1 " S'J ' (47) 

Taken at face value, eq. (|4.7| ) predicts infinite coupling constant at exponentially small mo- 
mentum scale, 

p = exp(— 27r/3A). 

This fact is a direct analog of the Landau pole in field theories without asymptotic freedom. 
In reality, it is a manifestation of the existence of a two-particle bound state in our model 
(recall that at least one bound state exists in every, arbitrarily weak, two-dimensional at- 
tractive potential). One can show that in our case, the energy of the bound state is of order 

Pi 



The 1 — > n amplitude is now 



n(n — 1) 

A n {t) = AT(l-^ 2 ■ (4.8) 

Let us note in passing, that we may compare our result for the case of large n, n ^> 1, to 
that obtained in ref. [^SJ for (2+1) dimensions by a different method. In the case when At 
is small, Xt <S 1, our formula in fact reproduces the result of ref. [^S 
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A n = A^ e exp(^-nh 



However, if At is comparable to one, eq. (|4.8|) does not coincide with that obtained in ref. 
H28fl . We consider this fact as a counter-argument to the claim of ref. 



C. O(N) theory 

The application of the technique described above to the multi-component if 4 model 
requires a slight modification. The lagrangian of the O(N) model, 

L = ^{d^a){d^ a ) - -^PofPa ~ j(<Pa<PaY ', (4.9) 

(a = 1 ... N is the internal index) contains one coupling constant A. However, if one tries 
to write down the most general O(N) symmetric effective non-relativistic lagrangian, one 
sees that there may exist two different effective couplings, g\ and g 2 , which correspond to 
the two possible structures of the potential term, 

L eff = ^ a id^ a - \{d&\){d& a ) - gi *t*l* b * b - 2g 2 ^ a ^ a %. (4.10) 

g\ is the low-energy elastic scattering scattering amplitude, singlet in the s-channel, while 
gi determines amplitudes, singlet in t- and w-channels. From the initial lagrangian ( ]4.9| ) 
one obtains the bare value of the coupling constants, 

A 

9i = 92= g • (4.11) 

The fact that the bare values of g\ and g 2 are equal to each other is the remnant of the 
crossing symmetry of our intial lagrangian (|4.9| ). However, the evolution equations for the 
running coupling constants are different, 

dgi{t) __ l rA7 2 



dt n 



-[Ngi{t)+A gi {t)g 2 {t)l (4.12) 
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In order to simplify the RG equations, let us introduce, instead of gi(t), a linear combination 
of the two couplings, 



9(t) = 9i(t) + -g 2 (t). 



The RG equation for g(t) is simpler than that of g\(t), eq. (|4.12|) , 



dg(t) N 2f 



dt n 



(4.14) 



The solution to the RG equations, eqs.( 4.13] ) and ( 4.14J ), which satisfies the initial condition 
( |4.11| ), can then be found, 



9 « = ( 1 + f)i( 1 + (7V+2) ^)' 



In analogy to the simple ip A case, the production of n soft bosons from an initial particle 
with isospin a can be described by an effective operator A n ^ a (W b W b p n ~ 1 '' 2 (n must be odd). 
After some calculations we obtain the following RG equation for A n (t), 



dA n (t) _ 1 
dt 7i 



^(N + n- l)g(t) + (n - l) 2 (l - 1 ) .</,(/) 



A n {t). 



(4.15) 



Having substituted the formulas for ^(^) an d d(t) into eq.( [4.15| ) and solved it, one obtains 
the dependence of the 1 — > n amplitude on the logarithm of the characteristic momentum 
of final particles, 



A n {t)=AT[l + {N + 2)^-t 



1 + S* 



(4.16) 



where, as before 



\p. 



2 Ve 



Recall that the regime which we are describing is A — ► 0, Xt ~ 1, provided other parameters 
such as n and iV are fixed. 
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Eqs.( |4.6|) , ( [4.8| ) and (|4.16|) are our final leading logarithmic expressions for the amplitudes 
1 — > n. As discussed in Introduction, they have the exponential form at Xt <C 1. 

To end up this section, let us consider the large N limit and show that the result of ref. 



25| can be reobtained in (2+1) dimensions by our technique. Consider the case when n is 
much smaller than both the number of boson flavors N and the inverse coupling constant. 
Then eq.( [4.16|) reduces to a simpler formula, 



Mt) = AT 



NX ' 

87T 



-(n-l)/2 



(4.17) 



Recalling that the coupling constant enters into the tree amplitude through the factor 
\(n~i)/2^ one can rewr jte eq. (|4.17|) in the following form, 



A n (t) ~ n\\%- l)/ \ (4.18) 

where 

is just the renormalized singlet scattering amplitude. One sees that for a small number of 
final particles, n <C N, leading order in 1/N result is given just by the tree-level formula 
where the coupling constant is replaced by the renormalized one. This is precisely the result 
of ref. |25| in the particular case of (2+1) dimensions. 

However, eq.( [4.18|) is valid only when n <C N. If the number of final particles is com- 
parable to N, the effect of loops is obviously not a simple renormalization of the coupling 
constant. One finds from eq.( [4.16|) that the correction to the large-iV result, eg. ( [4.171 ), 
is proportional to n 2 /N. When the number of final particles becomes comparable to the 
number of their spieces, the 1/N expansion becomes unreliable. One can expect that the 
breakdown of the 1/N expansion is not a peculiar feature of (2+1) dimensions but holds 
also in (3+1)- and higher-dimensional theories. 

So, we see that the renormalization group is a poweful mean for investigating the mul- 
tiparticle amplitudes in (2+1) dimensional scalar field theory at and around the threshold. 
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The exact formula for the amplitude, if ever be found, must incorporate the information 
obtained here by making use of the renormalization group equations. 

V. ni -► n 2 PROCESSES. 

A. Heuristic arguments. 

Let us consider the scattering of n\ virtual particles of equal energies ui\ into n 2 virtual 
particles of equal energies uj 2 in (p 4 theory. The amplitude of this process is given by the 
following path integral, 

A 

I D(pe iS ( f <p(x,t)e iuJlt dxdtY 1 (f <p(x,t)e- iu2t dxdtY^ 
Let us study the case when 

^1,2 

m, 2 = — 

where V\^ 2 are fixed as A — > 0. The amplitude fl5.1|) can be written in the equivalent form, 

A 
^■n 1 ^n 2 — 

(5.2) 

[Dtp exp (iS + n 1 \n(f <p(x, t)e iuJlt rfx dt) + n 2 \n(f <p(x, t)e~ iuJ2t rfx dt) J 
The change of variables 

transforms the integral (~ |5.2[) into an apparently saddle point form 

Ani ^ 12 = (5.3) 

A -(n 1+ n 2 )/2 j D( p exp i (iSty) + Vl i n Q 0( X; t )e iuJlt dx dt) v 2 In (/ <f>(x, t)e~ iuJ2t dx dt)) 

which indicates that the amplitude is exponential, 

A n ^ n2 oc A-(" 1+ " 2 )/ 2 exp(^(z/ 1 ,z/ 2 )) 
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Equivalently, this expression can be written in a way conjectured in Introduction, 

4>™ = yfc\yfitie* F( - Vl ' Va) (5.4) 

where 

F(z/i, z/ 2 ) = W(^i, i/ 2 ) - — m i/i - — In v 2 + — + — 

Similar observations can be made in the case when final and/or initial particles are on 
mass shell; a natural way to proceed in that case is to use the coherent state representation. 

Note that the above observations do not take into account possible cancellations in the 
pre-exponential factor, which, in fact, do appear, at least at the tree level p9|-pl[]. 



We have not been able to convert the above observations into a proof for the following 
reason. In general, eq.( |5.1|) contains disconnected graphs, so it does not, in fact, correspond 
to the connected amplitude. One may try to single out the connected amplitude by imposing 
a constraint that n\ and n 2 are coprime numbers (up to a common divisor two). However, 
in that case the conservation of energy, n\OJ\ = ^2^2? would mean that u\ and uo 2 entering 
eq. (|5.3|) are not arbitrary, but their ratio is a ratio of two large natural coprime numbers (the 



same is true for the ratio V\/v-z). In this way large numbers enter the exponent in eq.(pl 
implicitly, in addition to explicit 1/A. Thus, the saddle point nature of the integral ([5.3|) is 
questionable. 

Nevertheless, we expect the functional form of the amplitude, eq. ( |5.4|) , to be correct. We 
check this conjecture by numerical calculations of the tree amplitudes in the next subsection. 

B. Numerical calculation of tree amplitudes. 

To construct the algorithm for numerical calculations of the tree amplitudes we make 
use of the following observation made in ref . |JT | . 



Let us consider amplitudes of scattering of n\ virtual particles into n 2 real particles when 
all particles are at rest. To avoid disconnected graphs, we keep n\ and n 2 being coprime 
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numbers up to one common divisor two. Let us construct an iteration solution to source-free 
classical space-independent field equation, 

(p + ip + Xip 3 = 

(we again set the mass equal to one) with the first iteration 

0(0) = Zl + Z2} (5.5) 

z x = Cie^\ z 2 = ( 2 e- i " 2t , (5.6) 

uo 2 = 1 when we are interested in scattering into real particles. The iteration procedure is 
defined by the following equation 

^) +y ,W = _ A(v9 3 )(fc -l) (57) 

where (y? 3 )^ -1 ) is of (k — l)-th order in A and is expressed through (p(°\ . . . , ip^ k ~ l \ At the 
/-th step (/ = (rii + n 2 )/2 — 1), in addition to oscillating terms, a peculiar term appears for 
the first time, and the amplitude A ni ^ n2 is determined by its coefficient, 

(J) = te -it A 2 (UJ) I C n lC n 2 -l + osciUating terms (5. 8 ) 

n x ! (n 2 - 1)! 2 x 

This algorithm was adapted for computer calculations. When one expands cp® in powers 
of exponents z\ and z 2 , eq.( |5.7| ) enables one to obtain the coefficients of the expansion of 
ip^ algebraically. This induces an efficient numerical procedure for the calculation of the 
amplitude A ril ^ ri2 . 

Since we are calculating only tree amplitudes, we can determine the power of A in the 
resulting amplitudes directly, by counting the number of vertices in the diagrams, 

< e ln 2 oc A^- 1 (5.9) 

By comparing eqs.(p^) and (p.9|), one finds that the function F should be a homogeneous 
function of the first order, so it should have the following form, 
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F tree (Am, \n 2 ) = \( ni + n 2 )$(z/) 
where v = 112/111, and the entire tree amplitude can be rewritten as follows, 

A^tn, = K - !) ! (^2 - l)!A^- 1 e( ni+na >*< na / ni)+0 < n °> (5.10) 

So, to check our conjecture at the tree level, we calculate the function 

n( "- n - )=bg (n 1 -1S:-i) l (5 - n) 

and verify that in the large-n limit, this function Q is linear at ii\ at given v. So, we should 
calculate Q at different ii\ and at fixed v and then observe that at large n\ it tends to a 
linear function at each v. 

In real calculations, we cannot fix v exactly: if we choose a particular m, then at given 
v the "number of outgoing particles", vn\, would not be, in general, integer. So, for fixed u, 
as we vary m, we can only choose an integer n 2 in such a way that B2/V11 is close to v, and 
112 is coprime to n\\ the larger m, the closer 712/711 to v. To have better precision, for each 
rii we have used three values of 112 giving n 2 /ni closest to v (two from above and one from 
below) and then made the interpolation to the chosen value of v. The functions Q u (iii), 
obtained in this way, are shown in fig. [X0| (unbroken theory) and fig. [TT] (broken theory). 



Figs. [Ty and |Tl| are consistent with Q v (ni) being linear function of n\ at large n\ at 
various v. This means that the function $(711,712) entering eq.( |5.10|) indeed depends only 
on n 2 /iii, i.e., our conjecture (|5.4|) is valid at the tree level. 



VI. CONCLUSIONS 

In complete analogy to the instanton-like processes in high energy collisions, there 
emerges an intriguing situation with calculations of the amplitudes of processes with large 
number of outgoing particles (n = 0(1/A)) in scalar theories. The results of this paper 
strongly suggest that the corresponding cross section is of the exponential form, 

an(xe ±G(Xn,XE) (61) 
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where E is the total center-of-mass energy. This form indicates that the exponent G may 
be calculable in a semiclassical manner, but the actual calculational technique is presently 
lacking. The perturbative calculations can at best provide the evaluation of first terms in 
the expansion of G in (Xn). This expansion blows up at (Xn) ~ 1, so the perturbation 
theory is of no use for studying the most interesting values of n. 

In the instanton-like case, several proposals have been put forward which might enable 
one to calculate the exponent for the cross section. These include the Landau technique for 
the calculation of the semiclassical matrix elements Hl6j|l7| and the study of multiparticle 



initial states [|l2j -|T4||. If either of these approaches works for multiparticle production in 
scalar theories, the exponent in eq. (|6.1|) should be independent of a particular choice of the 
initial state (whether it contains one or several particles, whether these particles are virtual 
or real, etc.). We hope that this expectation can be checked by making use of the techniques 
developed in this paper. 

It might happen that the exponent G in eq. ( |6 . 1| ) is negative at all n and E, so that the 
multiparticle cross sections are always exponentially small. Even in that case the calculation 
of the exponent would be of substantial interest as it would require the development of novel 
methods of quantum field theory. 

The authors are indebted to M.B.Voloshin for pointing out the peculiar infrared prop- 
erties of scalar theories in (2+1) dimensions, for numerous useful discussions and help 
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FIGURES 




FIG. 1. Process "few — ► many" 





FIG. 2. Recurrence relation for tree amplitudes 
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FIG. 3. Feynman rules for calculating multiparticle amplitudes 
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FIG. 4. The one-loop graph (a) and a representation (a') of its leading singularity 
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FIG. 5. The two-loop graphs with their symmetry factors (a-f) and their leading singularities 
(a'-f) 
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FIG. 6. A disconnected graph giving subleading contribution at large n. 




FIG. 7. The graph contributing to the lowest-order beta function in the effective 
non-relativistic theory. 
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FIG. 8. One- and two-loop graphs that contribute to the 1 — > 3 amplitude in the infrared 




FIG. 9. A graph that makes contribution to RG equation for 1 — ► n amplitude. There are 
n(n — l)/2 diagrams of this type, corresponding to different possiblities to choose a pair of final 
particles to rescatter. 

FIG. 10. Numerical results for n\ — ► ri2 tree amplitudes in the unbroken (p 4 theory (see 
Postscript file) 

FIG. 11. Numerical results for n\ — ► ri2 tree amplitudes in the broken <p 4 theory (see Postscript 
file) 
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